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Abstract 

The current attempt is aimed to honor the first centennial of 
Johannes Diderik van der Waals (VDW) awarding Nobel Prize in 
Physics. The VDW theory of ordinary fluids is reviewed in the first 
part of the paper, where special effort is devoted to the equation of 
state and the law of corresponding states. In addition, a few mathe- 
matical features involving properties of cubic equations are discussed, 
for appreciating the intrinsic beauty of the VDW theory. A the- 
ory of astrophysical fluids is shortly reviewed in the second part of 
the paper, grounding on the tensor virial theorem for two-component 
systems, and an equation of state is formulated with a convenient 
choice of reduced variables. Additional effort is devoted to particu- 
lar choices of density profiles, namely a simple guidance case and two 
cases of astrophysical interest. The related macroisothermal curves 
are found to be qualitatively similar to VDW isothermal curves below 
the critical threshold and, for sufficiently steep density profiles, a crit- 
ical macroisothermal curve exists, with a single horisontal inflexion 
point. Under the working hypothesis of a phase transition (assumed 
to be gas-stars) for astrophysical fluids, similar to the vapour-liquid 
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phase transition in ordinary fluids, the location of gas clouds, stellar 
systems, galaxies, cluster of galaxies, on the plane scanned by reduced 
variables, is tentatively assigned. A brief discussion shows how van 
der Waals' two great discoveries, namely a gas equation of state where 
tidal interactions between molecules are taken into account, and the 
law of corresponding states, related to microcosmos, find a counter- 
part with regard to macrocosmos. In conclusion, after a century since 
the awarding of the Nobel Prize in Physics, van der Waals' ideas are 
still valid and helpful to day for a full understanding of the universe. 
keywords - galaxies: evolution - dark matter: haloes. 

1 Introduction 

One century ago (1910), the Nobel Prize in Physics was awarded to Johannes 
Diderik van der Waals (hereafter quoted as VDW). In his doctoral thesis 
(1873) the ideal gas equation of state was generalized for embracing both the 
gaseous and the liquid state, where these two states of aggregation not only 
merge into each other in a continuous manner, but are in fact of the same 
nature. With respect to ideal gases, the volume of the molecules and the 
intermolecular tidal forces were taken into account. 

The VDW equation was later reformulated in terms of reduced (dimen- 
sionless) variables (1880), which allows the description of all substances in 
terms of a single equation. In other words, the state of any substance, defined 
by the values of reduced volume, reduced pressure, and reduced temperature, 
is independent of the nature of the substance. This result is known as the 
law of corresponding states. 

The VDW equation of state in dimensional and reduced form, served 
as a guide during experiments which ultimately led to hydrogen (1898) and 
helium (1908) liquefaction. The Cryogenic Laboratory at Leiden had devel- 
oped under the influence of VDW's theories. For further details on VDW's 
biography refer to specialized textbooks (e.g., Nobel Lectures 1967). 

The current paper has been written in honor of the first centennial of 
VDW awarding Nobel Prize in Physics. The ideal and VDW equation of 
state, both in dimensional and reduced form, are reviewed, and a number 
of features are analysed in detail, in Section [2l Counterparts to ideal and 
VDW equations of state for astrophysical fluids, or macrogases, are briefly 
summarized and compared with the classical formulation in Section [31 The 
discussion and the conclusion are drawn in Section m 
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2 Equation of state of ordinary fluids 



Let ordinary fluids be conceived as fluids which can be investigated in lab- 
oratory. The simplest description is provided by the theory of ideal gas, 
where the following restrictive assumptions are made: (i) particles are iden- 
tical spheres; (ii) the number of particles is extremely large; (iii) the motion 
of particles is random; (iv) collisions between particles or with the wall of 
the box are perfectly elastic; (v) interactions between particles or with the 
wall of the box are null. 

The equation of state of ideal gases may be written under the form (e.g., 
Landau and Lifchitz, 1967, Chap. IV, §42, hereafter quoted as LL67): 



where p is the pressure, V the volume, T the temperature, the particle 
number, and k the Boltzmann constant. 

In getting a better description of ordinary fluids, the above assumption 
(v) is relaxed and tidal interactions between particles are taken into con- 
sideration. The VDW generalization of the equation of state of ideal gases, 
Eq. ([1]), reads (van der Waals, 1873): 



where A and B are constants which depend on the nature of the particles. 
More specifically, the presence of an attractive interaction between particles 
reduces both the force and the frequency of particle-wall collisions: the net 
effect is a reduction of the pressure, proportional to the square numerical 
density, expressed as AlN/V)"^. On the other hand, the whole volume of the 
box, V, is not accessible to particles, in that they are conceived as identical 
spheres: the free volume within the box is — NB, where B is the volume 
of a single sphere. For further details refer to specific textbooks (e.g., LL67, 
Chap. VII, §74). 

The isothermal (T = const) curves for ideal gases are hyperbolas with 
axes, p = ^V, conformly to Eq. ([T]). In VDW theory of real gases, the isother- 
mal curves exhibit two extremum points below a threshold, which reduce to 
a single horisontal inflexion point when a critical temperature is attained, 
as shown in Fig.lH Well above the critical isothermal curve, T ^ Tc, the 
trends exhibited by ideal and VDW gases look very similar. Below the criti- 
cal isothermal curve, T < T^, the behaviour of VDW gases is different with 
respect to ideal gases and, in addition, the related isothermal curves provide 
a wrong description within a specific region where saturated vapour and liq- 
uid phases coexist. Further details are shown in Fig. El Above the critical 



pV = kNT 



(1) 
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Figure 1: Isothermal curves related to ideal (left panel) and VDW (right 
panel) gases, respectively. Isothermal curves (from bottom to top) corre- 
spond to T/Tc = 20/23, 20/22, 20/21, 20/20, 20/19, 20/18. No extremum 
point exists above the critical isothermal curve, T/T^ = 1. 



isothermal curve (T = T^.) the trend is similar with respect to ideal gases. 
Below the critical isothermal curve and on the right of the dashed curve, 
the supersaturated vapour still behaves as an ideal gas. Below the critical 
isothermal curve and on the left of the dashed curve, the liquid shows little 
change in volume as the pressure rises. Within the bell-shaped area bounded 
by the dashed curve, the liquid phase is in equilibrium with the saturated 
vapour phase. A diminished volume implies smaller saturated vapour frac- 
tion and larger liquid fraction at constant pressure, and vice versa. The 
VDW equation of state is no longer valid in this region. The dashed curve 
(including the central branch) is the locus of intersections between VDW and 
real isothermal curves, the latter being related to constant pressure where 
liquid and vapour phases coexist. The dotted curve is the locus of VDW 
isothermal extremum points. 
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A specific [T/T^ = 20/23) VDW and corresponding real isothermal curve, 
are represented in Fig. [21 The VDW isothermal curve and the real isothermal 
curve coincide within the range, V < Va and V > Ve- The VDW isother- 
mal curve exhibits two extremum points: a minimum, B, and a maximum, 
D, while the real isothermal curve is flat, within the range, Va < V < Ve- 
Configurations related to the VDW isothermal curve within the range, Va < 
V < Vb (due to tension forces acting on the particles yielding superheated 
liquid), and Vd < V < Ve (due to the occurrence of undercooled vapour), 
may be obtained under special conditions, while configurations within the 
range, Vb < V < Vb, are always unstable. The volumes, Va and Ve, cor- 
respond to the maximum value in presence of the sole liquid phase and the 
minimum value in presence of the sole vapour phase, respectively. 

The surfaces, ABC and CDE, are equal, as first inferred by Maxwell (e.g., 
Rostagni, 1957, Chap. XII, §19). The VDW and real isothermal curves 
represented in Fig. [3] being related to the same temperature, T, the cycle, 
ABCDECA, is completely both isothermal and reversible, and the work, W, 
performed therein cannot be positive to avoid violation of the second law of 
the thermodynamics. The cycles, ABCA and CDEC, occurring in counter- 
clockwise and clockwise sense, respectively, are also completely both isother- 
mal and reversible. Accordingly, VFabcdeca = W/abca — W^cdec < 0. A sim- 
ilar procedure, related to the reversed cycle, ACEDCBA, yields PVacedcba = 
W^CEDC - M^CBAC < 0. Then PFabcdeca = W^acedcba = 0, which implies 
W^ABCA = WcDEC = W^CEDC = W^CBAC and, in turn, the equality between the 
related surfaces. For further details refer to specific textbooks (e.g., LL67, 
Chap. VIII, §85). 

In order to simplify both notation and calculations, it is convenient to 
deal with (dimensionless) reduced variables (e.g., Rostagni, 1957, Chap. XII, 
§16; LL67, Chap. VIII, §85). To this aim, the first step is the knowledge of the 
parameters related to the critical point, Vc,Pc,Tc. Using the VDW equation 
of state, Eq. ([2]), the pressure and its first and second partial derivatives, with 
respect to the volume, read: 

kNT 

A— ■ iV = const ; (3) 
V ^ 

kNT N^ 

+ 2^T7T ; (4) 




{V-NBy 1/3 
2kNT ^ N"^ 

6^771 ; (5) 



- NBf 1/4 

where the domain is 1/ > NB, V = NB is a vertical asymptote, and p = 
is a horisontal asymptote. The critical isothermal corresponds to the highest 
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temperature allowing a liquid phase, which occurs therein only at the crit- 
ical point. The critical isothermal curve exhibits neither a minimum nor a 
maximum, which are replaced by a horisontal inflexion point coinciding with 
the critical point. Accordingly, {dp/dV)v^,T^ = 0, {d'^v/dV'^)v^^T^ = 0, and 
Pc = kNTjlVc - NB) - AN^/V^. The solution of the related system is: 

K = 3iV5 ; (6) 
8 A 1 



1 A 

PcVc 

kNT, 



(9) 



c 



where, in general, the compressibility factor, Z = pV/{kNT), defines the 
degree of departure from the behaviour of ideal gases, for which Z = 1, 
according to Eq. ([T]). For further details refer to specific textbooks (e.g., 
Rostagni, 1957, Chap. XII, §20; LL67, Chap. VIII, §85). 
With regard to the reduced variables: 

r-v.^ ^-i- r-l- (1°' 

the ideal gas equation of state, Eq. ([1]), and the VDW equation of state, 
Eq. ([2]), reduce to: 

iI>Y=It ; (11) 

{^+l,]{f-\) = lr ; f>\: (12) 



V 3/ 3 ' ' 3 
and Eqs. ([3]), (jlj), and ([5]), reduce to: 

^~ 3f - 1 ~ f2 5 
( 2Af 6 



3 



drLr (3f-l)^ f 
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(13) 
(14) 

(15) 



where, for assigned ^, the domain of the function, ^()^), is Y > 1/3, Y = 
1/3 is a vertical asymptote, and ^ = is a horisontal asymptote. In the 
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special case of the critical point, )^ = 1, /' = l,j(=l, the partial derivatives 
are null, as expected. 

The extremum points, via Eq. (1141) . are defined by the relation: 

/(f) = ^^7^=^' (16) 

which is satisfied at the critical point, as expected. The function on the 
left-hand side of Eq. f lT6|) has two extremum points: a minimum at )^ = 1/3 
(outside the physical domain) and a maximum at )^ = 1, where ^ = 1. 
Accordingly, Eq. f lT6|l is never satisfied for ^ > 1, which implies no extremum 
point for related isothermal curves, as expected. The contrary holds for 
^ < 1, where it can be seen that the third-degree equation associated to 
Eq. ([Tl]) has three real solutions, related to extremum points. One lies outside 
the physical domain, which implies Y ^ 1/3- The remaining two are obtained 
as the intersections between the curve, f{f\ expressed by Eq. (fT6|) . and the 
straight line, y = ^, keeping in mind that /(1/3) = 0, /(I) = 1, and 
lim^^+oo/(f) = 0. 

The third-degree equation associated to Eq. ( TM|) . may be ordered as: 

f3-9af2 + 6af -a = ; (17a) 

« = 47 ' (l^b) 
with regard to the standard formulation (e.g., Spiegel, 1968, Chap. 9): 

+ a\X^ + a^x + as = ; (18) 
the discriminants of Eq. fll7ap are: 

Q = ^^^ = a(2-9a) ; (19) 

9aia2 - 27a3 - 2a? a(l - 18a + h^c?) 
^ = 54 = 2 ' (^°) 

D = Q^ + R^='^^^^ ; (21) 

where D = in the special case of the critical isothermal curve ( = 1, a = 
1/4), D < for < 1, and D > for ^ > 1. Accordingly, three coincident 
real solutions exist if D = 0, three (at least two) different real solutions if 
D < 0, one real (outside the physical domain) and two complex coniugate if 
D>0. 
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The three real solutions {D < 0) may be expressed as (e.g., Spiegel, 1968, 
Chap. 9): 



-(y cos I vr + - 



f2 = 2^-gCOS (^TT 

f3 = 2^/^cos(7r + - + ^ 



3 
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47r' 
3" 



3^^ 



6* = arctan 



R 



(22a) 
(22b) 
(22c) 
(22d) 



where ai = 
a = 1/4, g 
f c = fo = 



—9a and, in the special case of the critical isothermal curve, 
-1/16, D = 0, which implies fo = min(fi, fs, fg), f a = f b = 
fE = max(fi, fa, fs)- In the special case, ;r 0, Eq. (fTTaj) 
reduces to a second-degree equation whose solutions are fi = fo2 = 1/3, 
while the related function is otherwise divergent as a — t- +oo. In general, 
the extremum points of VDW isothermal curves ( ^ < 1) occur at f = fe 
(minimum) and f = f) (maximum), f^ < fj. As ^ — t- 0, fe — ?■ 1/3, 
Yd — ^ +00, where, in all cases, l/3< )^<1< f). 

The two areas defined by the intersection of a generic VDW isothermal 
curve ( < 1) and related real isothermal curves (see Fig.[3D, are expressed 
as: 



Wi 
W2 



nVc rVc 

/ pcdV - pdV = pcVc 



Va 
Vc 



pdV 



Va 
Ve 



Vc 



Pc dV = pcVc 



•Ye 



fA 



Yc 



Ya 



^df 



/ "^df -^c(fE- fc) 

JYc 



(23a) 
(23b) 



and the substitution of Eq. (IT^ into 
integrals. The result is: 



Wi 

PcVc 

PcVc 
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^c(fc-fA)-3ri.3y^_l 
3' 3fc- 



allows explicit expressions for the 

3fc-l , 3(fc- fA) 



3(fE- fc) 

fcfE 



fAfc 
^c(fE- fc) 



(24a) 
(24b) 



and the condition, Wi = W2, after some algebra reads (Caimmi 2010, here- 
after quoted as CIO): 



8 ^ 



In 



3fE-l 

3fE-fA'"3fA-l 

8 



fAfE 



(25) 



where, for a selected isothermal curve, the unknowns are i>c = i>A = i>Ei Yki 
and Ye- 

The reduced volumes, ^c, ^e, see Fig.|3l may be considered as inter- 
sections between a VDW isothermal curve ( ^ < 1) and a horisontal straight 
line, ]/) = 1^(7, in the (Oy^^) plane. In other words, ^a, Yc, Te, are the real 
solutions of the third-degree equation: 

f'-(l + l^]f' + ^f-^ = o ; (26) 



3 3 



which has been deduced from Eq. (1131) . particularized to ^ = ^c- The re- 
lated solutions may be calculated using Eqs. fl22l) . The last unknown, i^c, is 
determined from Eq. (1251) . 

An inspection of Fig.[3] shows that the points, A and E, are located on the 
left of the minimum, B, and on the right of the maximum, D, respectively. 
Keeping in mind the above results, the following inequality holds: < 
)^ < 1 < yt) < ^E, which implies further investigation on the special case, 
Yc = 1- The particularization of the VDW equation of state, Eq. (|T3l) . to the 
point, C = Ci, assuming Yci = 1, yields: 



^=^£i±l. (27) 



and Eq. (126|) reduces to: 



f3-(l + 26)f2 + 36f -6 = ; (28a) 

6=-^ ; (28b) 

with regard to the generic third-degree equation, Eq. flTK]) . the three solutions, 
Xi, X2, Xs, satisfy the relations (e.g., Spiegel, 1968, Chap. 9): 

xi + X2 + X3 = -ai ; (29a) 

X1X2 + X2XS + X3X1 = a2 ; (29b) 

X1X2X3 = -as ; (29c) 

where, in the case under discussion: 

ai = — 1 — 26 ; a2 = 36 ; 03 = — 6 ; (30a) 

xi = fA; X2 = fc, = l; xs=fE; (30b) 



and the substitution of Eqs. fl30l) into two among fl29|) yields: 



fA = 6-V623^ ; (31a) 
fE = b + VW^ ; (31b) 
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and the combination of Eqs. fl27|) . (I28bp . and (1^ produces: 



fA 



1-2^1- y 



T < 1 



(32a) 



1+ yi- y 

4;r-3 



(32b) 



which, together with Yci = 1; are the abscissae of the intersection points be- 
tween a selected VDW isothermal curve in the (0)^^) plane and the straight 
line, 1^ = ^Ci, in the special case under discussion. 

The substitution of Eqs. fl32|) into fl25|) . the last being related to the real 
isothermal curve, yields: 



which holds only for the critical isothermal curve, ^ = 1. Accordingly, the 
abscissa of the intersection point, C, between a selected VDW isothermal 
curve and related real isothermal curve, see Fig.[3l cannot occur at = 1 
unless the critical isothermal curve is considered. Then the third-degree equa- 
tion, Eq. fl26|) . must be solved in the general case by use of Eqs. fl22|) . The 
results are shown in Tab.dl where the following parameters (in reduced vari- 
ables) are listed for each VDW isothermal curve, see Fig. El the temperature, 
^, the lower volume limit, Y^^, for which the liquid and vapour phase coexist; 
the extremum point (minimum) volume, fs', the intermediate volume, fc, 
for which the pressure equals its counterpart related to the corresponding 
lower and upper volume limit, for which the liquid and vapour phase coexist; 
the extremum point (maximum) volume, Yd] the upper volume limit. Ye, for 
which the liquid and vapour phase coexist; the extremum point (minimum) 
pressure, ^s; the pressure, = = ^e, related to the horisontal real 
isothermal curve; the extremum point (maximum) pressure, The locus 
of the intersections between VDW and real isothermal curves is represented 
in Fig. [2] as a trifid curve, where the left, the right, and the middle branch 
correspond to Ya, Ye, and Yc, respectively. The common starting point 
coincides with the critical point. The locus of the VDW isothermal curve 
extremum points is represented in Fig.[2]as a dotted curve starting from the 
critical point, where the left and the right branch corresponds to minimum 
and maximum points, respectively. 

A fluid state can be represented in reduced variables as (Y, T), where 
one variable may be expressed as a function of the remaining two, by use 




(33) 
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Table 1: Values of parameters, ^, Tb, Tc, Td, Te, ^b, ^c, ^d, within 
the range, 0.85 < < 0.99, using a step. A/" = 0.01. Additional values are 
computed near the critical point, to increase the resolution. The true value of 
the reduced temperature on the last row is ;r = 0.9999 or 10 ;r = 9.999. All 
values equal unity at the critical point. Index captions: A, C, E - intersections 
between VDW and real isothermal curves; B - cxtrcmum point of minimum; 
D - extremum point of maximum. Extremum points are related to VDW 
isothermal curves, while their real counterparts are flat in presence of both 
liquid and vapour phase. To save aesthetics, 01 on head columns stands for 
unity. 
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of the reduced ideal gas equation of state, Eq. flTT]) . or the reduced VDW 
equation of state, Eq. f ll2p . The formulation in terms of reduced variables, 
Eqs. (fTOj) . makes the related equation of state universal i.e. it holds for any 
fluid. Similarly, the Lane-Emden equation expressed in polytropic (dimen- 
sionless) variables, describes the whole class of polytropic gas spheres with 
assigned polytropic index, in hydrostatic equilibrium (e.g., Chandrasekhar 
1939, Chap. IV, §4). 

The states of two fluids with equal ( f ■, i>-, T)i defined as corresponding 
states. The mere existence of an equation of state yields the following result. 

Law of corresponding states. Given two fluids, the equality between two 
among three reduced variables, Y> in^plies the equality between the 

remaining related reduced variables i.e. the two fluids are in corresponding 
states. 

The law was first formulated by van der Waals in 1880. For further details 
refer to specific textbooks (e.g., LL67, Chap. VIII, §85). 

3 Equation of state of astrophysical fluids 

Let macrogases be defined as two-component fluids which interact only gravi- 
tationally. For assigned density profiles, the virial theorem can be formulated 
for each subsystem, where the potential energy is the sum of the self potential 
energy of the component under consideration, and the tidal energy induced 
by the other one. The virial theorem for subsystem can be expressed as a 
macrogas equation of state in terms of dimensionless variables, Xy, Xp, Xt, 
related to axis ratio, mass ratio, virial (i.e. self + tidal) potential energy 
ratio, respectively. The result is (CIO): 

XpXvFx{Xp,Xv)=XT ; (34a) 
X„ = m^ ; Xv = - ; Xt = (P ; (34b) 

y 

where the function. Fx, depends on the selected density profiles, m is the 
(outer to inner component) mass ratio, y is the (outer to inner component) 
axis ratio along a generic direction, (p is the (outer to inner component) 
virial energy ratio, and the density profiles are restricted to be homeoidally 
striated. The variables, Xy, Xp, Xt, play a similar role as the volume, 
the pressure, and the temperature, for ordinary fluids. Accordingly, Xy, Xp, 
Xt, may be defined as macrovolume, macropressure, and macrotemperature, 
respectively. For further details refer to the parent paper (CIO). 
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Macroisothermal curves on the (OXyXp) plane exhibit a similar trend 
with respect to VDW isothermal curves on the (OVp) plane, with two main 
differences. First, no critical point occurs for sufficiently mild density profiles, 
where all macroisothermal curves are characterized by two extremum points, 
one maximum and one minimum. Second, a critical macroisothermal curve 
appears for sufficiently steep density profiles, above (instead of below) which 
macroisothermal curves exhibit extremum points. For further details refer 
to the parent paper (CIO) and an earlier attempt (Caimmi and Valentinuzzi 
2008). 

The last inconvenient may be avoided turning Eq. (1341) into the following: 

Y,YvFy{Y„Yv) = Yt ; 

FyiY„Yv) = FxiXj„Xv) 

as suggested in the parent paper (CIO). 

The existence of a phase transition moving along a selected macroisother- 
mal curve, where the path is a horisontal line ("real" macroisothermal curve) 
instead of a curve including the extremum points ( "actual" macroisothermal 
curve), must necessarily be assumed as a working hypothesis, due to the 
analogy between VDW isothermal curves and macroisothermal curves. Un- 
like the VDW equation of state, Eq. the theoretical macrogas equation of 
state, Eq. (I35ap . is not analytically integrable, which implies the procedure 
used for determining a selected macroisothermal curve, must be numerically 
performed. 

The main steps are (i) calculate the intersections, Yy^, Yy^, Yve? ^Va < 
Yvc < ^Ve5 between the generic horizontal line in the (OYyYp) plane, Yp =const, 
and the theoretical macrogas equation of state, within the range, l^g <Yp< 
Yp^, where B and D denote the extremum points of minimum and maximum, 
respectively; (ii) calculate the area of the regions, ABC and CDE; (iii) find the 
special value, Yp = Yp^, which makes the two areas equal; (iv) trace the real 
macroisothermal curve as a horisontal line connecting the points, {Yyj^, V^^), 
{Yy^,Yp^), (lVg,l^g), Ypj^ = Yp^ = = Yp^. For further details refer to an 
earlier attempt (CIO). 

The procedure related to point (ii) above is rather cumbersome and should 
be performed again with the new variables, Yy, Yp, and Y^, with respect to an 
earlier attempt (CIO). For this reason, the current paper shall be restricted 
to theoretical macroisothermal curves and related extremum points. In order 
to preserve the analogy with ideal and VDW gases, the tidal potential energy 
shall be excluded and included, respectively, in the formulation of the virial 



(35a) 

Yt = ^ ; (35b) 

(35c) 
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theorem and related equation of state. The following cases shall be dealt 
with: UU macrogases, where no critical point occurs; HH macrogases, where 
the critical point occurs; HN/NH macrogases, where the critical point occurs. 

In presence of the critical point, Eq. ( 135|) may be translated into reduced 
variables, as: 



X,YyFy{Xp.Xv)^^ = XT ; (36a) 

% = ^ ; rv^ = ^ ; rT = ^ ; (36b) 

^Pc -^K: 

Fy(yp,yv) = Fy(>;yp^,>V>Vj ; (36c) 

where Yp^, Yy^^ Yt^, are the values of the variables related to the critical point. 
The counterpart of Eq. fl36ap for ideal macrogases reads: 

Y^YyGy{rM^^ = rT ; (37) 

where Gy^Yp, Yv) is the expression of Fy{Yp, Yv) where the interaction terms 
are omitted. For further details refer to an earlier attempt (CIO). Ac- 
cordingly, the equation of state for ideal macrogases where Gy{Xp,Yv)Yp^/ 
{Yvc^Tc) = 3/8, coincides with its counterpart related to ideal gases, con- 
formly to Eq. (fTTl) . 

Macroisothermal curves related to lUU (tidal potential energy excluded) 
and AUU (tidal potential energy included) macrogases, are plotted in Fig.lH 
left and right panel, respectively, for values of the macrotemperature. It = 
20/23, 20/22, 20/21, 20/20, 20/19, 20/18, from bottom to top. The coordi- 
nates, Yy, Yp, Ft, may be conceived as normalized to their fictitious critical 
counterparts, Yy^ = 1, Yp^ = 1, Yt^ = 1 (CIO). The comparison with ideal 
and VDW gases, plotted in Fig.[Tl shows a similar trend, except the ab- 
sence of a critical macroisothermal curve, above which the extremum points 
disappear. 

Macroisothermal curves related to IHH (tidal potential energy excluded) 
and AHH (tidal potential energy included) macrogases, are plotted in Fig. [5], 
left and right panels, respectively, for infinitely extended subsystems and val- 
ues of the reduced macrotemperature, Yt = ^t/^t^ = 20/23, 20/22, 20/21, 
20/20, 20/19, 20/18, from bottom to top. The general case of bounded sub- 
systems makes only little changes. The comparison with ideal and VDW 
gases, plotted in Fig.[T|, shows a similar trend where macroisothermal curves 
are more extended along the horisontal direction with respect to isothermal 
curves. 
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Macroisothermal curves related to IHN/NH (tidal potential energy ex- 
cluded) and AHN/NH (tidal potential energy included) macrogases, are plot- 
ted in Fig. [HI left and right panels, respectively, for infinitely extended sub- 
systems and values of the reduced macrotemperature, Yt = Yt/Yt^ = 20/23, 
20/22, 20/21, 20/20, 20/19, 20/18, from bottom to top. The general case 
of bounded subsystems makes only little changes for AHN/NH macrogases, 
while the scale change tends to disappear for IHN/NH macrogases. The 
comparison with ideal and VDW gases, plotted in Fig.[Tl shows a similar 
trend where macroisothermal curves are more extended along the horisontal 
direction with respect to isothermal curves, and the occurrence of a scale 
difference for ideal macrogases. The last is due to a mass divergence for in- 
finitely extended N density profiles, which makes tidal effects higly increase. 

The comparison between the VDW critical isothermal curve and its coun- 
terparts related to HH and HN/NH macrogases is shown in Fig. [3 The bro- 
ken curve is the same as in Fig.O Accordingly, the vapour and the liquid 
phase of ordinary fluids coexist within the bell-shaped region bounded by 
the broken curve. Both HH and HN/NH macroisothermal curves are more 
extended along the horisontal direction with respect to VDW isothermal 
curves, which implies a more flattened counterpart of the above mentioned 
bell-shaped region. The critical point belongs to all curves. 

4 Discussion and conclusion 

Tidal interactions between neighbourhing bodies span across the whole ad- 
missible range of lengths in nature: from, say, atoms and molecules to galax- 
ies and clusters of galaxies i.e. from micro to macrocosmos. Ordinary fluids 
are coUisional, which makes the stress tensor be isotropic and the velocity 
distribution obey the Maxwell's law. Tidal interactions (electromagnetic 
in nature) therein act between colliding particles (e.g., LL67, Chap. VII, 
§74). Astrophysical fluids are collisionless, which makes the stress tensor 
be anisotropic and the velocity distribution no longer obey the Maxwell's 
law. Tidal interactions (gravitational in nature) therein act between a single 
particle and the system whole (e.g., CIO). 

In both cases, an equation of state can be formulated in reduced vari- 
ables: the VDW equation for ordinary fluids and an equation which depends 
on the density profiles for astrophysical fluids. For sufficiently mild density 
profiles, macroisothermal curves are characterized by the occurrence of two 
extremum points, similarly to isothermal curves where a transition from liq- 
uid to gaseous phase takes place, or vice versa. For sufficiently steep density 
profiles, a critical macroisothermal curve exhibits a single horisontal inflex- 
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ion point, which defines the critical point. Macroisothermal curves below 
and above the critical one, show two or no extremum point, respectively, in 
complete analogy with VDW isothermal curves. In any case, the existence 
of an equation of state in reduced variables implies the validity of the law of 
corresponding states for macrogases with assigned density profiles. 

For astrophysical fiuids, the existence of a phase transition must neces- 
sarily be assumed as a working hypothesis by analogy with ordinary fiuids. 
The phase transition has to be conceived between gas and stars, and the 
{OYvYp) plane may be divided into three parts, namely (i) a region bounded 
by the critical macroisothermal curve on the left of the critical point, and 
the locus of onset of phase transition on the right of the critical point, where 
only gas exists; (ii) a region bounded by the critical macroisothermal curve 
on the left of the critical point, the locus of onset of phase transition on the 
left of the critical point, and the vertical axis, where only stars exist; (iii) a 
region bounded by the locus of onset of phase transition, and the horisontal 
axis, where gas and stars coexist. The locus of onset of phase transition, 
not shown in Fig. [7] for reasons explained above, is similar to its counterpart 
related to ordinary fiuids, represented by the bell-shaped curve in Fig. [71 but 
more extended along the horisontal direction. 

In this view, elliptical and SO galaxies lie on (ii) region unless hosting hot 
interstellar gas, and the same holds for globular clusters; spiral, irregular, and 
dwarf spheroidal galaxies lie on (iii) region, and the same holds for cluster 
of galaxies; gas clouds in absence of star formation lie on (i) region, and the 
same holds for hypothetic galaxies with no stars. 

In conclusion, van der Waals' two great discoveries, more specifically a 
gas equation of state where tidal interactions between molecules are taken 
into account, and the law of corresponding states, related to microcosmos, 
find a counterpart with regard to macrocosmos. After a century since the 
awarding of the Nobel Prize in Physics, van der Waals' ideas are still valid 
and helpful to day for a full understanding of the universe. 
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Figure 2: Same as in Fig.[T] (right panel), where the occurrence (within the 
bell-shaped area bounded by the dashed curve) of saturated vapour is consid- 
ered. Above the critical isothermal curve (T = T^) the trend is similar with 
respect to ideal gases. Below the critical isothermal curve and on the right 
of the dashed curve, the gas still behaves as an ideal gas. Below the crit- 
ical isothermal curve and on the left of the dashed curve, the liquid shows 
little change in volume as the pressure rises. Within the bell-shaped area 
bounded by the dashed curve, the hquid phase is in equilibrium with the sat- 
urated vapour phase. A diminished volume implies smaller saturated vapour 
fraction and larger liquid fraction at constant pressure, and vice versa. The 
VDW equation of state is no longer valid in this region. The dashed curve 
(including the central branch) is the locus of intersection between VDW and 
real isothermal curves, the latter being related to constant pressure where 
liquid and vapour phases coexist. The dotted curve is the locus of VDW 
isothermal extremum points. 
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Figure 3: A specific (T/Tc = 20/23) VDW and corresponding real isothermal 
curve. The above mentioned curves coincide within the range, V < Va 
and V > Ve- The VDW isothermal curve exhibits two extremum points: a 
minimum, B, and a maximum, D, while the real isothermal curve is flat within 
the range, Va ^ V < Ve- Configurations related to the VDW isothermal 
curve within the range, Va < V < Vb (due to tension forces acting on 
the particles yielding superheated liquid), and Vb < V < Ve (due to the 
occurrence of undercooled vapour), may be obtained under special conditions, 
while configurations within the range, Vb < V < V^,, are always unstable. 
The volumes, Va and Ve, correspond to the maximum value in presence of 
the sole liquid phase and the minimum value in presence of the sole vapour 
phase, respectively. The regions, ABC and CDE, have equal area. For further 
details refer to the text. 
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Figure 4: Macroisothermal curves related to lUU (left panel) and AUU (right 
panel) macrogascs. respectively. Macroisothermal curves (from bottom to 
top) correspond to Ft = 20/23, 20/22, 20/21, 20/20, 20/19, 20/18. No 
critical macroisothermal curve exists, above which the extremum points dis- 
appear. The coordinates, Yy, Yp, Y^, may be conceived as normalized to 
their fictitious critical counterparts, Yy^ — 1, Yp^ — 1, Yt^ — 1. 
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Figure 5: Macroisothcrmal curves (^p = Y^jYp^ vs. "j/y = ly/lV^) related 
to IHH (left panels) and AHH (right panels) macrogases, respectively, for 
infinitely extended subsystems. Macroisothermal curves (from bottom to 
top) correspond to Yt = ^t/^t, =20/23, 20/22, 20/21, 20/20, 20/19, 20/18. 
The general case of bounded subsystems makes only httle changes. 
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Figure 6: Macroisothermal curves = Yp/Yp^ vs. Yy = Yy/YyJ 

related to IHN/NH (left panels, to be noted the scale difference) and 
AHN/NH (right panels) macrogases, respectively, for infinitely extended 
subsystems. Macroisothermal curves (from bottom to top) correspond to 
= Yt/Yt, =23/20, 22/20, 21/20, 20/20, 19/20, 18/20. The general case 
of bounded subsystems makes only little changes for AHN/NH macrogases, 
while the scale difference tends to disappear for IHN/NH macrogases. 
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Figure 7: Comparison between VDW critical isothermal curve (full), HH 
critical macroisothermal curve (dotted) and HN/NH critical macroisothermal 
curve (dot-dashed). With regard to ordinary fluids, the vapour and the hq- 
uid phase coexist within the bell-shaped region bounded by the dashed curve 
and, in addition, ly = V^^p = P- More extended (along the horisontal direc- 
tion) bell-shaped regions are expected for HH and HN/NH macroisothermal 
curves. The critical point belongs to all curves. Different letters denote the 
expected location of different astrophysical systems. Caption: EG - eUiptical 
galaxies; SO - lenticular galaxies; SG - spiral galaxies including barred; IG 
- irregular galaxies; DS - dwarf spheroidal galaxies; GC - globular clusters; 
CG - clusters of galaxies; WC - wholly gaseous clouds i.e. in absence of star 
formation; WG - (hypothetical) wholly gaseous galaxies i.e. in absence of 
star formation. 
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